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Abstract. Small, but significant non-random fluctuations in event-by-event average 
Pt have been observed in Au+Au collisions at ^/sjvat = 200 GeV by the PHENIX 
collaboration. These are consistent with being caused by correlations due to jets at 
large pr, where the measured suppression must be included to reproduce the centrality 
dependence of the non-random fluctuations. 



1. The event-by-event average pr distribution is not a Gaussian, it's a 
Gamma distribution. 

The single particle inclusive px distribution averaged over all particles in all events 
in a p-p experiment (inclusive) or in a given centrality class in an A+A experiment 
(semi-inclusive) is usually written in the form: 

_^_ = 6V^ or ^ = b 2 p T e- h ^ . (1) 
PrdpT dpx 

Equation [T] represents a Gamma distribution with p = 2, where, (px) = p/b, a PT / (p T ) = 

an d typically 6 = 6 (GeV/c)^ 1 for p-p collisions. The 'inverse slope parameter' 

T = 1/b is sometimes referred to as the 'Temperature parameter'. 

For events with n detected charged particles with magnitudes of transverse 

momenta, pr v the event-by-event average px, denoted M PT is defined as: 

1 « i 

M PT =p^=-Y J PT 1 = -E Tc . (2) 
n ~{ n 

For the case of statistical independent emission, where the fluctuations are purely 
random, an analytical formula for the distribution in M PT can be obtained 
assuming negative binomial (NBD) distributed event-by-event multiplicity, with Gamma 
distributed semi-inclusive pt spectra. pQ The formula depends on the 4 semi-inclusive 
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parameters (n), 1/k, b and p which are derived from the means and standard deviations 
of the semi-inclusive pr and multiplicity distributions, (n), a n , (pr), o~ PT : 

f(y)= J2 /nbd(w, 1/k, (n)) f T (y,np,nb) , (3) 

71 ^min 

where y = M PT . For fixed n, and purely random fluctuations, the mean and standard 
deviation of M PT follow the expected behavior, (M PT ) = (pr), o~m vt = o~ PT /y/n. In 
PHENIX, equation El is used to confirm the randomness of mixed-events (Figure |IJ). 




Figure 1. PHENIX mixed-event distribution for the 0-5% centrality class (data 
points) compared to equation [3] (curve). 

2. Measurement of non-random fluctuations in PHENIX 

Mixed-events are used to define the baseline for random fluctuations of M PT in 
PHENIX (21 E] • This has the advantage of effectively removing any residual detector- 
dependent effects. The event-by-event average distributions are very sensitive to the 
number of tracks in the event (denoted n or N trac k s ), so the mixed event sample is 
produced with the identical N tra cks distribution as the data. Additionally, no two tracks 
from the same data event are placed in the same mixed event in order to remove any 
intra-event correlations in p?. Finally, (M PT ) must exactly match the semi- inclusive 
(pr)- As noted above, the randomness of M PT for the mixed-event sample is tested 
by comparison to equation El Figure Q shows the excellent agreement between the 
calculation and the mixed event M PT distributions for the 0-5% centrality class. The 
standard deviations, &m Pt , differ by less than 0.04%. This represents the maximum 
error from any effects introduced by the event mixing procedure. 

The measured M PT distributions for the data in two centrality classes for ^snn = 
200 GeV Au+Au collisions in PHENIX are shown in Figure 121 (data points) compared 
to the mixed-event distributions (histograms). The non-Gaussian, Gamma distribution 
shape of the M PT distributions is evident. The difference between the data and the 
mixed-event random baseline distributions is barely visible to the naked eye. The non- 
random fluctuation is quantified by the percent difference of the standard deviations of 
M PT for the data and the mixed-event (random) sample: 

„ &M PT ,data 0~M VT ,mixed , . 

Pp T = - . (4J 

&M PT , mixed 
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Figure 2. M PT for 30-35% and 0-5% centrality classes: data (points) mixed-events 
(histogram). 



The results are shown as a function of centrality represented by N part in Figure El 
compared to the previous PHENIX measurement at ^snn = 130 GeV. j2j The errors 
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Figure 3. F PT in % :(left) Au+Au, p+p, ^/IWn = 200 GeV, (right) Au+Au 130 GeV 



shown are systematic errors due to time-dependent detector variations. Comparatively, 
statistical errors are negligible. The systematic error is calculated from the rms variation 
of F PT from 10 independent subsets of the data. The improvement over the ^/snn = 130 
GeV data is due to 3 times larger solid angle (larger N tracks ), better tracking and more 
statistics. |3] 

The dependence of F PT on N part is striking. To further understand this dependence 
and the source of these non-random fluctuations, F PT was measured over a varying 
p T range, 0.2 GeV/c < p T < p™ ax (Figure EJ), where p™ ax = 2.0 GeV/c for the N part 
dependence. The increase of F PT with p™ ax suggests elliptic flow or jet origin. This was 
investigated using a Monte Carlo simulation of correlations due to elliptic flow and jets 
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Figure 4. F PT vs centrality and p™ ax compared to simulations. 



in the PHENIX acceptance. The flow was significant only in the lowest centrality bin 
and negligible (F PT < 0.1%) at higher centralities. Jets were simulated by embedding [at 
a uniform rate per generated particle, S pro b(N part )} p-p hard-scattering events from the 
PYTHIA event generator into simulated Au+Au events assembled at random according 
to the measured N tracks and semi-inclusive px distributions. This changed (px) and a PT 
by less than 0.1% . S pro b(N part ) was either constant for all centrality classes, or scaled 
by the measured hard-scattering suppression factor RAA{N par t) for pt > 4.5 GeV/c. |lj 
A value F PT = 2.06% for p-p collisions was extracted from pure PYTHIA events in the 
PHENIX acceptance in agreement with the p-p measurement. The value of S pro b(N part ) 
was chosen so that the simulation with S pro b(N part ) x RAA{N part ) agreed with the data 
at Np ar t = 182. The centrality and p™ 2 ™ dependences of the measured F PT match the 
simulation very well, but only when the Raa scaling is included. 

A less experiment-dependent method to compare non-random fluctuations is to 
assume that the entire F PT is due to temperature fluctuations of the initial state, with 
rms variation a T /(T): 012] 

«n)-l) o* T /(T)* _ P 4 
FpT ~ 2 olJ{p T )*-2 [{n) V> 2 • () 

This yields o~t/ {T)=1.8% for central collisions and 3.7% at the peak of F PT , which puts 
severely small limits on the critical-fluctuations that were expected. 
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